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Abstract. The origin of the prominent 'knee' in the cosmic ray energy spectrum at 
an energy of several PeV is still uncertain. A recent mechanism has shown promise, 
however; this involves particles from a very young pulsar interacting with the radiation 
field from a very young supernova remnant. The ensuing nuclear reaction of the 
particles with the photons by way of e = e~ production then causes the characteristic 
knee. 

In an earlier paper we argued that the mechanism would imply only one source of a 
very rare type - if it were to explain the spectral shape. Here we examine the mechanism 
in more detail and conclude that for even a single source to work its characteristics 
would need to be so unusual that the mechanism would not be possible for any known 
type of pulsar-supernova combination. 
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1. Introduction 

There are only two features of the cosmic ray energy spectrum which depart from 
smoothness: the 'knee' at a few PeV and the 'ankle' at a few EeV. The origin of both 
is still subject to argument. Here we examine the first feature. 

The 'knee', a comparatively sudden change of slope - a steepening - has been known 
since the pioneering work of Kulikov and Khristianssen [1] . Most workers have attributed 
it to the rather sudden reduction in Galactic trapping efficiency, the trapping being 
confinement by the Galactic magnetic field. The different atomic nuclei will have 'cut- 
offs' at Z-dependent energies and, when summed, the result will be a rather smooth 
(not sharp) change in the spectral slope; we term this the Galactic Modulation, GM, 
model. A minority of workers have gone along with our idea of the knee being due to 
a single, recent, local supernova remnant (SNR) ((2] and later publications) This is the 
single source model (SSM). Yet another, recent, model has been put forward by Hu et 
al. in Ref. [3] (to be referred to as I) in which particles from a rapidly rotating pulsar 
interact with radiation from the parent SNR. We refer to this model as PSNR. 

The aim of the present work is to follow up our earlier work ([!]), (to be referred 
to as II) in which we examined the PSNR Model as originally proposed in I to see if 
the knee caused by the particle -SNR photon interactions was a universal property of 
all pulsar, SNR pairs. We found that it was not but that there was hope for a single 
rare pair having the required properties. This would be identified as the single source 
of the SSM although the mechanism would differ from that in the original SSM which 
involved SNR alone. Here we examine the needed properties of the pulsar, SNR pair in 
more detail. 

2. The basis of the PSNR model. 

2.1. The characteristics of the interactions. 

As remarked already, the philosophy of the source process is a composite pulsar-SNR 
system with the pulsar-accelerated particles interacting with the photon field of the 
SNR. Under certain circumstances a sharp knee will result due to particle, photon — > 
e + e~ interactions. We start by examining the general properties of the interaction in 
terms of the mean energy of the photons. 

An important characteristic of the 'knee' in the energy spectrum is its 'sharpness', 
S. We follow Ref. j2] and define the sharpness as 

(log IE 3 ) 

d(\ogE) 2 ' 1 ' 

Figure [T] sets the scene for the attenuation length, R, of protons and iron interacting 
with black body photons. The energy density of the photons is 1 eV cm -3 for all 
but the Cosmic Microwave Background (CMB), where it is the usual 0.24 eV cm -3 . 
The minimum range is inversely proportional to the mean photon energy, T. (1 eV 
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corresponds to 1.16 x 10 4 K, ie a mean wavelength of 12,400 A). We note that the 
middle of the visible range, yellow light, has wavelength ~ 5000 A ie a mean energy of 
2.48 eV. 

The effect of the photon field on the energy spectrum of the emerging particles 
depends on the energy density of the photons at the time of their passage through 
the field, for a low energy density there is no effect and for a high energy density the 
interactions are catastrophic, with no particles emerging. As an illustration, with the 
object of finding S- values in the required range we choose a combination of energy 
density and 'trapping time' such as to give the minimum R in Fig. [1] for T = 1 eV 
(R = 3 x 10 30 cm). This must be equal to c/ E(t) t dt, where E(t) is the energy density 
at time t. The limits, t and tx, are discussed in the next section. So as to maximise 
the PSNR 'effect' we choose the smallest value of to an d the largest t\ permissible. 

2.2. The limiting time parameters for the interactions and the two models. 

The lower limit to the interaction time, t , is set by energy losses by the pulsar- 
accelerated particles on the environmental gas density. This latter is enhanced very 
considerably by the ejecta from the SN; it is this material, too, that obscures the optical 
radiation for a period of days. For a typical massive progenitor of a Type II SN (and a 
'magnetar') the mass ejected is of order 5 M & [5] The expansion velocity at short times 
after the expansion can be found by reference to model calculations only; a commonly 
quoted value is ~ c/10 Some support for a value of this order comes from [5A], who 
measured c/30 some 30 years after the expansion of a shell remnant (43.31, +592) in 
Messier 82. 

If the density of the ejecta were constant, an attenuation length for protons of 120 g 
cm -2 would be achieved in ~ 1.6 x 10 6 s, neglecting the ISM swept up. If concentrated 
in a thin shell the time needed would be ~ 10 6 s. In fact, there will be considerable 
dumpiness, reducing this time by at least an order of magnitude; indeed, observations 
of (much older) SNR show very variable densities. Here, as a rather extreme limit (a 
factor 30 smaller) we adopt t = 3 x 10 4 s. A measure of support comes from the fact 
that, occasionally amongst the very many SNR light curves one encounters 'light break 
through' at such early times. 

Turning to the maximum time, t±, its' derivation is easier. Inspection of Figured] 
shows that for a photon 'temperature' of leV and an energy density of 1 eV cm -3 , 
interactions will cease when R ~ 10 30 cm. This can be converted to an appropriate radius 
for the environment of a SNR of known total photon energy and its time dependence. The 
highest energy output achievable is ~ 10 43 erg s _1 [6j yielding in 10 6 s : 10 49 erg. 

We now come to an important problem: the extent to which the pulsar-accelerated 
particles are trapped in the remnant before escaping. Alternatives are adopted here: 

(i) the particles are trapped by the tangled magnetic fields and diffuse out with the 
remnant (Model 1) and 
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(ii) the particles are not trapped but stream out from the pulsar with the speed of light 
(Model 2) 

Predictions will be made for each case and then common problems will be addressed. 

3. Results for Model 1. 

3.1. Model 1 - particle diffusion. 

SNR commonly expand at a velocity of ~ c/10 as already remarked and an expansion 
time of 10 6 s yields a radius of 3 x 10 15 cm. Assuming that the CR are trapped in the 
remnant for this time the required total distance traveled for the photon energy achieved 
corresponds to R ~ 10 30 cm for 1 eV cm -3 . The energy density for our remnant follows 
as ~ 5 x 10 13 eV cm" 3 . 

To summarise: our (rather extreme) limits are to = 3 x 10 4 s and t\ = 10 6 s; the 
corresponding SNR radii are 10 14 cm and 3 x 10 15 cm. 

3.2. Model 1 - The derived spectrum after the interactions. 

Calculations have been made using the above parameters, for protons in the first 
instance. The results are shown in Figure [21 for both the spectral reduction factors, 
/, and sharpness values, S. In order to study the sensitivity to the mean photon energy, 
the calculations are made for both T = 1 eV and 2 eV. 

It should be noted that the S- values are upper limits in that we ignore the particles 
which are accelerated after the interactions with the SNR-photons have effectively 
ceased. Such particles will dilute the knee and reduce S. 

Figure [3] gives the corresponding plots for Helium; it will be seen that the S- values 
are rather similar - the main difference being, as expected, in the energy at which S is 
a maximum. 

3.3. Model 1 - The ensuing knee position 

Inspection of Figures [2] and [3] shows that the energies for S m3X are, for protons: E*(p): 
0.31 PeV (1 eV) and 0.15 PeV (2 eV) and for helium nuclei: E*(H.e): 1.2 PeV (1 eV) 
and 0.62 PeV (2 eV) 

The ratio of the energies for S max for helium and protons are seen to be ~ 4, as 
would be expected since it is energy per nucleus that is important in determining the 
threshold energy for photon, particle interactions. 

Comparison can now be made with experiment. The energy of the knee in the CR 
energy spectrum is model-dependant and values range from about 3-8 PeV. There is 
evidence that Helium nuclei predominate at the knee [7], [Hj). If so, there is clearly a 
factor 10 difference in the S^ax-energies, between observed and expected, for T = 1 eV. 
A value nearer T = 0.1 eV (wavelength ~ 120, 000 A: for infra-red) is required. This is 
a serious drawback for Model I. 
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3.4- Model 1 - A variant to fit the CR spectral knee. 

Figure H] shows the situation where T is reduced to 0.13 eV so that the proton knee 
(peak) is now at E*(p) ~ 2 PeV. The Helium knee (peak) is at E*(He) ~ 7.8 PeV. Both 
can be regarded as, perhaps, just acceptable. 

Two aspects need consideration, now. The first concerns the position of the iron 
knee (peak), at 100 PeV. This is close to the 70 PeV where recent work [8] gives rather 
strong evidence for a peak from an analysis of the results from the EAS arrays. This is 
a good feature of the model. The second relates to the need for T ~ 0.13 eV; with the 
conversion 1 eV corresponds to 1.16 x 10 4 K, the effective temperature is in the region 
of 1,500 K, ie in the infrared, as already remarked. Although there is, inevitably, an 
infrared flux, from visible radiation absorbed by the SNR material, its intensity is too 
low to give sufficient energy density for our purpose. 

3. 5. Model 1 - The effect of a non-uniform SNR temperature. 

Inevitably, the temperature of the remnant will fall with time and a reduction of S will 
occur. Figure [5] shows the situation for the likely case of a 'conventional SNR' (eg [9]) 
where the temperature falls by a factor 2 over the operative period of particle-photon 
interactions. Interestingly, although finite, the reduction in S is not large. However, 
for more dramatic SNR, such as are required here, the change in temperature can be 
much bigger. In [6] temperatures versus time are given for the following progenitors: 
red supergiants, blue supergiants and Wolf-Rayet stars. The ratios of T(3 x 10 4 s) to 
T(10 6 s) are, respectively, (for typical stars within these types): 7.1, 8.4 and 8.4. The 
exponents of the power laws are -0.56, -0.61 and -0.61. It is true that the luminosity falls 
with time, too, but with much smaller exponents: -0.17, -0.35 and -0.35. In consequence, 
the reduction in S in practice for SNR from such progenitors will be much smaller than 
shown in figure [5j 

For young magnetars, potential SN for the mechanism considered here, the light 
curves are a rather strong function of the magnetic field of the magnetar and often 
continue for hundreds of days [8]. Predictions are not available for the short period 
from SN birth (0.3 to 10 days) needed here so that conclusions cannot be drawn. 

In summary, falling temperatures appear not to be a problem for the model. 

4. Results for Model 2 

Inspection of many published SNR 'luminosity curves' shows a rather wide spread but 
Figure shows the general form,for the case where the absorption by the SNR shell 
material is only 10%. Figure |6b shows the light curve for an absorption of 50%. The 
parameter for the lines is the 'cooling time', i.e. the mean exponential decay time of the 
radiation. 

Calculations have been made by us for a variety of values for the absorption 
efficiency, the cooling time and the time over which the source is active (the 'source 
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activity time'). Demanding that the sharpness should exceed S=2 gives the results shown 
in Figure [7J It will be noted that there is a minimum cooling time of 40 days and the 
source activity time must be greater than about 30+ days. 

5. Discussion. 

5.1. Parameters to be discussed. 

A number of parameters need examination, principally: 

the masses needed for injection, 
the total photon energy, 

the products of the photon-nucleus interactions and, most importantly, 
the characteristics needed for the pulsar. 

These parameters will be considered, in turn. 

5.2. The masses needed for injection 

It is apparent that the masses of nuclei needed: protons, Helium, Iron (?) are not 
expected as products of pulsar acceleration. Rather, proton-rich or iron-rich neutron- 
star surface elements would be expected. Although not an impossible situation, the 
needed composition does not add confidence to the model. 

5. 3. The total photon energy 

The total photon energy in the photon field adopted (which we assume effectively ceases 
after 10 6 s) is ~ 10 49 erg. To this should be added the later energy - which is the energy 
'seen' after the obscuration has ceased, but here we ignore this component which should 
not be large. The value can be compared with the total energy released in a typical 
SN: about 10 51 erg, of which perhaps 3 x 10 49 erg appears in visible light [9]. Thus, the 
energetics are acceptable. 

5-4- The products of the nucleus-photon interactions. 

If the whole of the CR spectrum were due to the PSNR mechanism then the overall 
change of slope of the primary CR spectrum would correspond to a transfer of energy 
from CR particles to gamma rays, via the electron pairs generated in the particle-SNR 
photon interaction. Inspection of the CR spectrum shows that the energy density of 
Galactic gamma rays would be about 2 x 10 -5 eV cm -3 , a value comparable with that 
detected in the diffuse gamma radiation, the latter usually being attributable to CR 
particles interacting with gas in the ISM. Such a situation would be 'dangerous' for the 
model but probably not catastrophic. In the more likely situation (in our view) that 
it is only the Single Source that produces the knee then copious gamma ray fluxes will 
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have been produced soon after the interaction (some few 10 5 y ago?). However, it is 
likely that these gamma rays will have 'been and gone' by now. 

The conclusion is that the non-observation of the gamma rays associated with the 
P-SNR interactions does not pose a threat to the model. 

5.5. The characteristics needed for the pulsars. 

The first characteristic is that the pulsar should be able to accelerate particles to 
rigidities of some lOPeV. Ref. [10] gives the expression 

Rmax = 6.6 X 10 12 B 12 P V 

Where B12 is the magnetic field in 10 12 Gauss, P is the Period in seconds and V is the 
rigidity in Volts. If we need H max = 10 16 V then, for B 12 = 3 (the mean, from 'model 
A' of [H]) a period of 45 ms is required. Interestingly, this is not too far from the birth 
period in the same model, which has P = 23 ms. 

The problem here is the likely age of pulsars of such characteristics. Using the work 
of [12] again, this is given by 

T = -^§- x 10 9 = 7 x 10 10 s 
2 B 12 

Such a time (1000 y) would be far too long to satify the requirements of Figure [3] 
(where the 'cooling time' is unlikely to be longer than ~ 100 d. 

Pulsars of much shorter period are required, for example, P ~ 1 ms and, say, 
B 12 = 3, for which T ~ 500 days. 

'Normal' pulsars (Lorimer, Model A, [TT]) have a standard deviation about the 
mean of cr(log-B) ~ 0.3 and will have a similar standard deviation for period [13], thus 
they will not satisfy the bill, even to the extent of 3 standard deviations from the mean. 

Magnetars [H] would, at first sight, appear to satisfy the requirements. They have 
extremely high magnetic fields, typically 6 x 10 14 G but, when measured, their periods 
are of order seconds. In [6] a theory is developed for the evolution of magnetars, in which 
the important birth period is in the range 10-740 ms, with a mean of 161ms. However, 
the mean field is given by (log I?) ~ 11.9. Clearly, the T-value will be, again, too long. 
Another factor militating against a magnetar as being the 'local' pulsar responsible 
for the knee (and the electron results) is the fact that their birth rate is quoted as 
~ 10~ 3 y _1 ([E]) furthermore none is in our vicinity; the catalogue, [33], gives the 
nearest as being some 3.5kpc away. In fact the most likely nearby SNR for the single 
source is a conventional one: MONOGEM [5]. 

6. Conclusions. 

That the suggested PSNR mechanism is very unlikely to work for even a single, rare 
'object' is provided by the following facts. 
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(i) The losses of energy by the pulsar-accelerated particles in collision with the very 
early SNR ejecta is expected to be very great. In our calculations we assumed 
that the ejecta was clumpy enough to allow such passage but this appears highly 
unlikely. 

(ii) For both models considered here the 'source activity time' is far too long for 
any known type of pulsar to be responsible; CR particles emitted after the SNR 
illumination has ceased will dilute the spectral sharpness to an unacceptable degree. 
This is intrinsic in Model 1 and, for Model 2, can be seen directly in Figure [TJfor a 
typical cooling time of 100 d the pulsar must be active for less than a year. 

(hi) Even if an ultra-rare pulsar of the required characteristics was proposed, the 
probability of its being close enough (and its contemporary presence undetected) 
is remote in the extreme. 

(iv) The mass composition of the ambient CR in general and the Simgle Source in 
particular are quite different from what would be expected from a pulsar. 
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Figure 1. Absorption length vs energy for protons and iron nuclei interacting with 
photon fields of mean 'temperature' as shown. The energy density of all but the CMB 
is 1 eV cm~ 3 




Figure 2. (a) The reduction factor, to be applied to the proton spectrum due 
to interactions with the photon field chosen (see text) for two 'temperatures', (b) 
Sharpness values derived from the above. 
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Figure 3. (a),(b) As for Figure 2 but for Helium nuclei. 




Figure 4. (a) /-values for the nuclei shown as a function of particle energy for T = 0.13 
eV, the temperature needed to give knee positions for protons and He near those found 
experimentally, (b) sharpness values derived from the above. 
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Figure 5. (a) /-values for a model in which the T- value fell smoothly from 2 eV to 
1 eV during the interaction period (full line) , in comparison with constant values of 2 
eV and 1 eV. (b)Sharpness values derived from the above. 




Figure 6. Luminisity for 10% (a) and 50% (b) absorption. 
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Figure 7. Model-2 in which the pulsar-accelerated particles travel in straight lines. 
The percentages refer to the fraction of the SNR radiation absorbed by the ejecta shell. 
'Source activity' relates to the time for which the pulsar is accelerating particles of the 
needed energy. The results are for helium. 



